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Abstract 

Charged tau leptons emerging in a long baseline experiment with a muon 
storage ring and a far-away detector will positively establish neutrino oscillations. 
We study the conversion of (Vn) and of I7e (z^e) to Ur or Vr for neutrinos from 
a 20 GeV muon storage ring, within the strong mixing scheme and on the basis 
of the squared mass differences which are compatible with all reported neutrino 
anomalies, including the LSND data. In contrast to other solutions which ignore 
the Los Alamos anomaly, we find charged tau production rates which should 
be measurable in a realistic set up. As a consequence, determining the complete 
mass spectrum of neutrinos as well as all three mixing angles seems within reach. 
Matter effects are discussed thoroughly but are found to be small in this situation. 
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1 Introduction 



Neutrino beams from a muon storage ring provide an ideal tool for the next round 
of experiments, aiming at establishing quantitatively oscillations and, possibly, CP- 
violation in the leptonic sector [Q, 0. The main reason for this is that both composition 
and spectrum of such beams are perfectly known and, in addition, the beam energy 
may be tuned. For instance, the neutrino beams from a storage ring such as the ones 
described in [|], |^ would contain an equal number of muon neutrinos and electron 
antineutrinos, or, depending on the sign of the parent muon, an equal number of muon 
antineutrinos and electron neutrinos, without any contamination from other neutrino 
species. 

The physics potential of a storage ring and a far-away detector, with regard to 
oscillations, was studied, e.g., in |P, §], and, with regard to CP- violation, in |^. In 
this paper we address the question of appearance of tau neutrinos in the processes 

z/^ — > Vr and — ^ , 

due to oscillations. We calculate the number of produced r+ and r~, taking into 
account matter effects on the neutrino beam on its way from the source to the detector, 
within the strong mixing scheme of three flavours only that we had proposed earlier 
10]. The latter assumption is receiving growing support by the fact that new data from 
KAMIOKANDE seem to disfavour the existence of a fourth, sterile, neutrino, while 
still showing the characteristic modulation of events as a function of zenith angle. As a 
consequence of the strong mixing which helps to populate strongly the Ur and V^- final 
channels from both the Vf, (z/g) and the (z7^) initial states, we find production 
rates one to two orders of magnitude larger than within a model ignoring the LSND 
data. 

In the next section we summarize some relevant formulae and collect the parameters 
extracted earlier from the analysis of all neutrino anomalies. We then describe and 
discuss matter effects in the Earth's crust and calculate the event rates for the 
appearance channels. The article ends with our results as well as some conclusions. 



2 Formulae and parameters 

In a scenario with three families of leptons, the mixing between three neutrino species 
is described by a conventional Cabibbo-Kobayashi-Maskawa matrix U relating flavour 
states to mass eigenstates. If neutrinos are Dirac particles, U has the form 



—iS " 

C12C13 512013 Si3e 

U = I — S12C23 — Ci2S23Si3e*'^ C12C23 — 512523^136*'^ ^23613 
S12S23 - Ci2C23Si3e*'^ -Ci2S23e*'^ - S12C23S13 C23C13 



1 



where Cij = cos% and Sij = sin 6ij. In the case they are of Majorana type, extra 
phases appear but their effects on neutrino observables are of order niy/ Ey and are 
generally negligible. 

The oscillation probabilities in vacuum take the form |0 

P{Vi Vj) = PcP-even{Vi Vj) + PcP-oddi^i Vj) (2) 
P{Vi Vj) = PcP-even{Vi ^ Vj) - PcP-oddi^i ^ Vj), (3) 

where 

Pcp-even{T^i ^ T^j) = Sij - ARbJH siu^ A12 - 4Re J^g siu^ A23 - 4ReJ^J sin^ A31, 

(4) 

Pcp-odd{vi Vj) = -SaijJ sin A12 sin A23 sin A31, 
with J the Jarlskog invariant and 

Jli^U,kUlU,f,Ul, A,,^AmjL/4^, cr,,^^:^,,, (5) 

k 

In order to account for all reported neutrino anomalies one needs only two squared 
mass differences. The LSND result can be understood in terms the mass difference 

AM^ := ml-ml^ 0.3 eV^ (6) 

The second mass difference is tuned so as explain the observed deficit of electron 
neutrinos coming from the sun, 

10"^ eV^ < Am^ < 10'^ (7) 

with 7712 = ^fn^i while the atmospheric neutrino data depend on both differences 
and are found in agreement with the prediction. Given the squared mass differences, 
the mixing angles are easily found to be (solution I) 

^12^35.5°, 023-27.3°, ^^13^13.1°. (8) 

This solution is favoured by the existing oscillation data and, thus, implies simultaneous 
and strong mixing of all three flavours. 
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3 Matter effects 



Of all neutrino species, only electron neutrinos can scatter elastically in the forward 
direction off electrons in matter, via charged current interactions. When the electron 
neutrinos oscillate into either muon or tau neutrinos, this introduces an additional 
term in the diagonal element of the neutrino flavour evolution matrix corresponding 
to Ve [|10[- is useful to define an effective mass term which stems from elastic 



scattering due to charged weak currents, 

a = 2V2GFn,E = 7.7 ■ 10" W ( (^) (9) 

ygr/cm'^y VGeVy 

where He is the electron number density in matter of density p and E is the neutrino 
energy. 

As is well-known matter effects become important only when a is comparable 
to, or larger than, the quantity A„2 = mj — m'j for some mass difference and neutrino 
energy. Given our neutrino mass spectrum and taking into account that for the Earth's 
crust p ~ 3gr/cm^, we are far from being in a range where matter effects would be 
dominant and could not be neglected. This would be the case, for example, in the case 
of — > z/e where large CP asymmetries are expected but will be masked by matter 
effects. In the case that we discuss here, matter effects are not expected to play a 
significant role but we will include them anyway. (In the case of electron neutrinos 
oscillating into tau neutrinos, and with our parameters, they represent at most a 2 % 
effect.) 

In order to exhibit the essential mechanisms we assume for a moment the two 
lightest neutrinos to be degenerate in mass. Indeed, this limiting case is close to the 
realistic situation where AM^ ^ Am^, cf. eqs. (P), (|^). In this case, the transition 
probabilities in vacuum for the "terrestrial" experiments depend only on three variables, 
i.e., 6*23, 6'i3 and AM^, as follows, 

= cos2(^23) sin2(2^i3) sin^ (^^) ' (10) 



P{u^^u,) ^ AUi,Ui, sin' 



= cos^(^i3) sin2(2^23) sin^ • (11) 

Note that in this limit of setting Am^ equal to zero the probabilities are independent 
of the angle 6*12. 
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When matter effects are included the above formulae are still valid provided one 
makes the replacements, 



with S23 unchanged. Here we have assumed that AM^ 3> a, Am^, a hierarchy which 
is respected by the mass spectrum we chose. 

From these formulae it is clear that the probability for a muon neutrino to oscillate 
into a tau neutrino in matter will not be different from its vacuum value. The reason 
for this is clear: the interpretation of the modulation of neutrino events with zenith 
angle reported by SuperKamiokande in terms of a muon neutrino oscillating into a 
tau neutrino either requires a nearly maximal 6'23, in two- neutrino mixing schemes, 
or still a sizeable one, in three-neutrino mixing schemes. (As noted above, the data 
strongly disfavour schemes with three active and one sterile neutrinos.) Furthermore, 
reactor experiments set a strict upper bound on 6'i3, giving sin^ 6*13 < .005. Thus, the 
factor cos'^(6'i3) in eq. (p^D, to a good approximation, is equal to 1. Likewise the factor 
sin^(2^23) in the same equation is also close to 1. Therefore, — > oscillations are 
favoured and the influence of matter effects on them is negligible. 

Introducing the terms containing Am^ does not modify this simple picture because 
they are suppressed by the huge gap between the two squared mass differences. It is 
important to notice that although we have assumed degeneracy between the two lightest 
neutrinos, this degeneracy is broken by matter effects which introduce an "effective" 
mass difference of the order 0(5^3 — l)/2. 

In the case of electron neutrinos oscillating into tau neutrinos, however, the story 
is different. In this case, the factor sin^(2^i3) in eq. (10) is a suppressing one and 
could compensate, at least partially, the gap in the squared mass differences. With 
our parameters, the contributions proportional to Am^ do not really compete with the 
one proportional to AM^ but they account for a sizeable correction. In our results 
presented below we use the full expression for the transition probabilities. 

At this point it is important to stress, that unlike our scheme, in schemes where 
only two anomalies are taken into account, disregarding the Los Alamos result, the 
typical mass differences are AM^ ^ 10^^ eV^ and Am^ ^ 10~^° eV^, or 10"^ eV^ to 
10~^ eV^, depending on whether in the vacuum solution the small-angle MSW solution 
or the large-angle MSW solution to the solar neutrino problem is chosen. Under these 
assumptions, AM^ is just in the appropriate range to exhibit sizeable matter effects 
so that these cannot be neglected. 

The conclusion so far is clear: If MiniBooNe confirms the LSND evidence for os- 
cillations, then AM^ is too large to cause a significant modification of the oscillation 




(13) 



(12) 
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probabilities due to matter effects. As we will see below, the prospects of discover- 
ing both z/^ — > Vj. and Vf. oscillations then are promising indeed. However, if 
KamLAND obtains a positive result in disappearance of electron neutrinos, corrobo- 
rating the large-angle MSW solution for solar oscillations, or if KARMEN definitely 
and conclusively excludes the LSND result, then AM^ will be in a range to produce 
a significant modification of the oscillation probability due to matter. In this case, 
the rates for tau appearance being lower by more than an order of magnitude, an ex- 
periment would be more difficult and would probably require a more intense neutrino 
source than that assumed here. 



= Ar^+ N^T ^^7/ / 9v, e P (^m ^ ^r) dE (15) 



4 Event rates 

The calculation of event rates for tau lepton (anti tau lepton) production from electron 
neutrinos (anti muon neutrinos) as a result of oscillations is straightforward. The total 
number of events in the two channels is given by 

lO'^ Na /-^^m 

n,- = NkT , / 9u. e <^ P {u, ^ u^) dE (14) 

10^ Na El fE,. 

where A'"^+ is the number of positive muon decays, a^^ is the charged current cross 
section per nucleon , P(z/e — > Ut) and P {u^ — > Ur) are the oscillation probabilities for 
neutrinos traveling inside the Earth taking into account matter effects. Nkx is the 
size of the detector in kilotons, 10^ Na is the number of nucleons in a kiloton, E^j, is 
the energy of the muons in the ring and = 5 GeV is a lower cut on the neutrino 
energies that we assume, for the sake of an example. For our numerical calculation, 
the energy spectrum (normalized to 1) of the neutrinos is taken to be 

g,^ = 12x\l-x) (16) 
g,^ = 2x\3-2x) (17) 

where x = Ey/ E^ is the fractional neutrino energy. It is straightforward to obtain the 
corresponding expression in the case of either electron anti neutrinos or muon neutrinos. 

Regarding the cross section, our calculation is based on the renormalization group 
improved parton model, focusing on the inclusive process //^(t^r) + N — ^ r~(r+) + 
anything. Note that in this case, unlike the charged current cross section for the muon, 
the terms proportional to the charged lepton mass are not negligible. The differential 
cross section is [O, 



4^ - {FM-y-^)^F.xy^±xFJy-y- 
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(18) 

ME\ ^\AE 2xJ 2^47] 

where we use the Bjorken scahng variables x = Q"^ /2Mv and y = v/Ey. Here —Q^ 
is the invariant momentum transfer between the incident neutrino and outgoing tau, 
u = Ey — Et is the energy loss in the laboratory frame while M and Mw are the 
nucleon and W-gauge boson masses respectively. 

The Callan-Gross relation 2xFi = F2 simplifies the above equation further. It is 
then dependent on only two form factors, F2 and F3, which are given in terms of the 
parton distributions by 

F2 = T.^(1^ + '^i) (19) 
F3 = E(9^-^0 (20) 



For our calculations we have used the MRST99 |T3| parton distribution. Although 
there does not seem to be a general consensus about the best way of combining the 
quasi-elastic resonance and the deep inelastic scattering cross sections at fixed energy 
to form the total cross section, a number of different approaches were proposed in the 
literature or were used in practice We adopt the simplest of them, and include 



only the deep inelastic part. In any case, the error one makes in the prediction due to 
uncertainties in the neutrino spectrum and by assuming a constant matter distribution 
in the Earth's crust is larger than the uncertainties in the cross section itself. 



5 Results and conclusions 

We consider a 20 GeV muon storage ring at CERN with a neutrino beam from one 
of its straight sections that points to the Gran Sasso underground laboratory, at a 
distance of 732 km. This is about the same distance as the one between FNAL and 
the Soudan mine. 

In order to mimic a nearly realistic experimental situation we set a lower cutoff in 
energy at 5 GeV and we assume a detection efficiency of 30%. Also, for two cases, 
production of r'^ from the process 17^ —>■ Vj. and production of from the process 
z/g — > i^^, we show our results in bins of 3 GeV so as to get a feeling for the number of 
events to be expected with our parameters. 

Fig. 1 shows the predicted spectrum for tau lepton appearance coming from v^. ^ fr-, 
while Fig. 2 shows the predicted spectrum of anti taus coming from V^j^ ^Vr- For the 
sake of comparison in Figs. 3 and 4 we plot the same observables for the choice of 
parameters reported in Ref which disregards the LSND data (note the different 
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scales in the ordinates!). Results similar to these were also reported in [|I^. Fig. 5, 
finally, shows our result for r~ and r+ production, grouped in bins of 3 GeV each. 

Provided the difference AM^ is in the range required by the LSND data, it is clear 
from Figs. 1, 2, and 5 that an experiment using a 20 GeV muon storage ring and a 
baseline of about 730 km between source and detector (corresponding to the distance 
CERN-Gran Sasso, or FNAL to the Soudan mine) would have a fair chance to see 
the appearance of charged tau leptons. The average oscillation probability could be 
measured with a statistical precision better than 3% and AM^ could be determined 
with a precision of about 3%. Note that since matter effects are small in this case, 
there is no additional uncertainty on this mass difference arising from an incomplete 
knowledge of the oscillation mode. 

All in all, with the choice of parameters that we obtained by a simultaneous ex- 
planation of all neutrino anomalies, a 10 kT detector some 732 km downstream would 
probe cos^(^23) sin^ (26*13) to very low values. When combined with measurements of 
t^fi oscillations (see for example |]T^) and with further results from atmospheric 



neutrinos, a precise determination of the complete mass spectrum of neutrinos and of 
all mixing angles seems possible. 
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Figure 1: Number of tau leptons detected in one year's time in a 732 km baseline 
and assuming 30% detecting efficiency. The solid line correspond to setting the CP 
violating phase 5 — while the dashed line correspond to 5 — 7r/2 
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Figure 2: Number of antitau leptons detected in one year's time in a 732 km baseline 
and assuming 30% detecting efficiency. The solid and the dashed curves (defined as 
in Fig.l) are indistinguishable. 
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Figure 3: Number of tau leptons detected in one year's time in a 732 km baseline 
and assuming 30% detecting efficiency for the parameters of ref. 0. Dashed and 
sohd hues as before. 
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Figure 4: Number of tau leptons detected in one year's time in a 732 km baseline 
and assuming 30% detecting efficiency for the parameters of ref. 0. The solid and 
the dashed curves are indistinguishable. 
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Figure 5: Number of tau (above) and antitau (below) leptons detected in one year's 
time in a 732 km baseline and assuming 30% detecting efficiency grouped in bins of 
3 GeV. 
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